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We explore whether an ultrasoft fermionic mode exists at extremely high temperature in Yukawa theory with massless 
fermion (coupling constant is g). We find that the fermion propagator has a pole at to = ±|p|/3 — i( for ultrasoft 
momentum p, where £ ~ g Ting , and the residue is Z ~ g . It is shown that one needs to take into account the 
asymptotic masses and the damping rate of hard particles to get a sensible result for such an ultrasoft fermionic mode; 
possible vertex correction turns out unnecessary for the scalar coupling in contrast to the gauge coupling. 



1. Introduction 

It is of basic importance to clarify the nature of quasi- 
particlc or collective excitations for understanding a many- 
r- — ipody system, in particular, in the low-energy regime. In 
'-^ a fermion-boson system such as quantum electrodynamics 
II (QED), quantum chromodynamics (QCD), and Yukawa 
Q-rtheory at so high temperature T that the masses of the 
Q particles are negligible, the average inter-particle distance 
i "is proportional to l/T and some collective effects can be 
expected in the soft momentum scale 1/gT with g being 
the coupling constant. Indeed, soft bosonic modes in the 
longitudinal as well as the transverse channels exist and 
are known as plasmon [1] , while the fermionic counter part 
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is known as plasmino [2, 3], both of which have masses of 
order gT. In the lower energy region, there exist hydro- 
dynamic modes of bosonic nature, which are actually the 
aero modes associated with the conservation of energy- 
momentum and charges. The energy hierarchy of bosonic 
and fermionic modes at high temperature may be summa- 
rized schematically as shown in Fig. 1. From the figure, 
one may have an intriguing but natural question whether 
such ultrasoft (< g 2 T) or zero modes can exist also in the 
fermion sector, possibly when the fermion system has a 
peculiar symmetry such as a chiral symmetry 

In this article, we argue and demonstrate that such a 
fermion mode can exist in this infrared energy region. Here 
we should mention that there have been some suggestive 
works for supporting the existence of such an ultrasoft 
fcrmionic mode at finite temperature. It was shown in a 
one-loop calculation [4, 5] that when a fermion is coupled 
with a massive boson with mass m, the spectral function 
of the fermion gets to have a novel peak in the far-low- 
cnergy region in addition to the normal fermion and the 
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plasmino, when T ~ m, irrespective of the type of boson; 
it means that the spectral function of the fermion has a 
three-peak structure in this temperature region. Very re- 
cently, the present authors [6] have suggested that such a 
three-peak structure may persist even at the high temper- 
ature limit in the sense m/T —> 0, for the massive vector 
boson on the basis of a gauge-invariant formalism, again, 
at the one-loop order. Thus, one may expect that the 
novel excitation may exist in the far-infrared region also 
for a fermion coupled with a massless boson, although the 
one-loop analysis admittedly may not be applicable at the 
ultrasoft momentum region. 

In supersymmetric models, the ultrasoft modes are dis- 
cussed as Nambu-Goldstonc fcrmions called phonino as- 
sociated with spontaneously breaking of supersymmctry 
at finite temperature, which is shown by using Ward- 
Takahashi identity and a diagrammatic technique [7, 8]. 
The discussion is extended to QCD at weak coupling [9], 
in which a supersymmetry is assigned at the vanishing 
coupling. The supersymmetry is explicitly broken by in- 
teractions; thus the phonino is not the exact zero mode 
but a pseudo-phonino. 

In this article, we analyze the fermion propagator in 
the ultrasoft energy region in Yukawa theory, which is 
a nonsupersymmetric model, because it is the simplest 
fermion-boson system and hence a good theoretical lab- 
oratory. We discuss the existence and the properties of 
the fcrmionic mode in the ultrasoft region using a similar 
diagrammatic technique in Rcfs. [8, 9, 10]. We deal with a 
massless fermion. In this case, the important point is that 
the fermion is chiral; the real part of the retarded fermion 
self-energy, E fl (w,0), at zero spatial momentum and the 
vanishing chemical potential, is an odd function of uj: 

Re£ fl (-o;,0) = -Re£ fl (w,0). (1) 

If the Re£ fl (u;,0) is a smooth function at uj = 0, then, 
Re E r (-oj, 0) = at u = 0, which implies Re G' 1 (w, 0) = 
at uj = 0; the spectrum has a peak at the origin provided 
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Figure 1: Schematic illustration of the energy hierarchy of bosonic 
and fermionic modes in a massless fermion-boson system at finite 
temperature: The vertical axis is the energy scale. The left and right 
part display the cases of bosonic and fermionic mode, respectively. 
The bosonic modes may be of scalar-, gauge-bosonic and hydrody- 
namic ones. No established fermionic mode had been known in the 
ultrasoft region. 
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Figure 2: Diagrammatic representation of the fermion self-energy 
in Eq. (2) at one-loop order. The solid and dashed line are the 
propagator of the fermion and the scalar boson, respectively. In 
Eq. (13), the fermion and the boson propagator in the internal lines 
are replaced by the dressed propagators given in Eq. (9). 



that the imaginary part of the fermion self-energy is not 
too large. This argument suggests that the existence of 
the ultrasoft pole may be a universal phenomenon at high 
temperature in the theory composed of massless fermions 
coupled with a boson, which includes not only Yukawa 
theory but also QED and QCD. 

It is, however, not a simple task to establish that 
fermionic modes exist in the ultrasoft region on a gen- 
eral ground beyond the one-loop order because of the in- 
frared divergence called pinch singularity [9, 10, 11, 12] 
that breaks a naive perturbation theory, as will be briefly 
reviewed in the next section. We remark that the same 
difficulty arises in the calculation of transport coefficients 
[10, 11] and the gluon self-energy [12]. 

We shall show that the retarded fermion propagator 
has a pole at to = ±|p|/3 — i( (£ ~ g 4 T\ng~ 1 is introduced 
later) with the residue Z ~ g 2 for ultrasoft momentum p. 

2. Ultrasoft Fermionic Mode 

As stated in the Introduction, we deal with a mass- 
less fermion (denoted by ip) coupled with a scalar field <f> 



through the interaction Lagrangian Ci = gipipfi a t high 
temperature T. We do not take into account the self- 
coupling of the scalar fields for simplicity. In this sec- 
tion, we calculate the fermion self-energy to obtain the 
fermion retarded Green function with an ultrasoft momen- 
tum p < g 2 T at high temperature. We first see that the 
naive perturbation theory breaks down in this case. Then, 
we shall show that a use of a dressed propagator gives a 
sensible result in the perturbation theory and that the re- 
sulting fermion propagator has a new pole in the ultrasoft 
region. 

The retarded self-energy in the one-loop level is given 

by 
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where Dq' s '(— fc) and Gq' S (p + k) are the bare propagators 
of the fermion and the scalar boson defined as 
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+ n B (fc )K(2^)sgn(fc°)<5(fc 2 ). (6) 
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Here, n F (k°) = l/(cxp(fc°/ T ) + 1) and n B (k°) = 
l/(exp(fc°/T) — 1) are the Fermi-Dirac and Bose-Einstein 
distribution functions, respectively. In the present analy- 
sis, we have preferentially employed the real-time formal- 
ism in Keldysh basis [14]. The diagrammatic represen- 
tation of Eq. (2) is shown in Fig. 2. Inserting Eqs. (3) 
through (6) into (2), we obtain 



S bare(p) = 



d A k 



n B (k°) 



X p2+ 2 p .r + ^ + ^ (2 - )Sgn(fc ° )J(fc2) 
x (27r)sgn(fc° + p°)S((k + p) 2 )], 

(7) 

where we have used the on-shell conditions for the bare 
particles, k 2 — 0, and (k + p) 2 = in Gq (p + k) and 
Dq (— k). Then, for small p, the self-energy is reduced to 
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where K{k) = (27r)sgn(fc°)(5(fc 2 )(n F (fc ) + n B (k )), which 
is independent of p. This approximation is equivalent to 
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the HTL approximation [2]. The HTL approximation is, 
however, only valid for p ~ gT, and not applicable in the 
ultrasoft momentum region. In fact, the retarded self- 
energy in the one-loop level obtained with use of the bare 
propagators is found to diverge when p — > 0, since the 
integrand contains 1/p-k, which singularity is the so called 
"pinch singularity" [9, 10, 11, 12]. 

The origin of this singularity is traced back to the use 
of the bare propagators because the singularity is caused 
by the fact that the dispersion relations of the fermion 
and the boson are the same and the damping rates are 
zero in that propagators. For this reason, one may sus- 
pect that this singularity can be removed by adopting the 
dressed propagators taking into account the asymptotic 
masses and decay width of the quasiparticles, as will be 
shown shortly. 

Noting that the leading contribution comes from the 
hard (k ~ T) internal and almost on-shell (k 2 w 0) mo- 
mentum 1 , we are lead to employ the following dressed 
propagators for the fermion and boson: 
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where m 2 



g 2 T 2 /8 and m 2 = g 2 T 2 /6 are the so called 
asymptotic masses of the fermion and the boson at A: 2 ~ 0, 
respectively [16, 17, 18]. The damping rates of the hard 
particles, Q and £&, are of order g 4 Tlng -1 . The logarith- 
mic enhancement for the damping rate of the scalar par- 
ticle is caused by the soft-fermion exchange, which is the 
similar enhancement to that of the hard photon [19, 20]. 
Using these dressed propagators, we have 
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for small p, where 5m 2 
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C/ + Cb- Here we have used the modified on-shell condition 
of the quasi-particles, k 2 — m 2 + 2iQk° = and k 2 — m\ + 
2i£fb = 0, to obtain the denominator of the integrand 
in Eq. (13). We neglected m b , m/, Cb, and Q in K(k), 
since the leading contribution comes from hard momenta 
k ~ T '. It is worth emphasizing that thanks to Sm 2 and 
C, S (p) given in Eq. (13) does not diverge in the infrared 
limit, p — > 0. 



1 Hcrc wc justify neglecting the case that the internal momenta are 
soft (k ~ gT) or smaller; for the soft momentum, the HTL resummed 
propagators [15] should be used. In these propagators, the dispersion 
relations of the fermion and the boson are different, so the pinch 
singularity does not appear. 



Here, we expand Eq. (13) in terms of 2(p-k+ik°() / Sr 
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with Z = 8(5?n 2 ) 2 /(g 2 n 2 T 4 ) and v = 1/3. Thus, we get 
the fermion propagator in the ultrasoft region as 
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Here we have decomposed the fermion propagator into the 
fermion number +1 and —1 sectors in the second line. 
These two sectors are symmetric under p <-» — p and 
v <-> —v, so we analyze only the fermion number +1 sector 
in the following: From Eq. (15), we find a new pole at 



uj = —v\p\ + i(. 



(16) 



Notice that the real part of u> is negative for the fermion 
sector, which suggests that this peak has an antifcrmion- 
hole-likc character like the antiplasmino [3]. The disper- 
sion relation of the real part, Rcw = — f|p|, is shown in 
Fig. 4 together with the HTL results [2, 3] for comparison, 
where the coupling constant is chosen as g = 0.1. The 
imaginary part of the pole reads 



(~g 4 T\ng-\ 



(17) 



which is much smaller than those of the normal fermion 
and the antiplasmino [16]. Since the real part and the 
imaginary part of the pole are finite for |p| ^ 0, this mode 
is a damped oscillation mode. The residue at the pole is 
evaluated to be 



Z = 



72tt 2 
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(18) 



which means that the mode has only a weak strength in 
comparison with those of the normal fermion and the an- 
tiplasmino, whose residues are order of unity. It is worth 
mentioning that such a smallness of the residue is actually 
compatible with the results in the HTL approximation: 
The sum of the residues of the normal fermion and the 
anti-plasmino modes obtained in the HTL approximation 
is unity and thus the sum rule of the spectral function of 
the fermion is satisfied in the leading order. Therefore, 
one could have anticipated that the residue of the ultra- 
soft mode can not be the order of unity but should be of 
higher order. 

Equations (15) through (18) for Yukawa theory are ob- 
tained for the first time and constitute the main results of 
the present work. We emphasize again that the usual HTL 
approximation [2, 3] is inapplicable in this ultrasoft energy 
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Figure 3: The spectral function in the fermion sector, Eq. (19), as a 
function of energy ui at zero momentum. The coupling constant is 
set to g = 0.1. 



Table 1: The respective coupling orders of the three fcrmionic modes. 
The momentum of the ultrasoft mode is set to <C g 2 T, while the 
momenta of the HTL results [1, 16] are set to zero. 
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region, and hence such an ultrasoft peak in the spectral 
function of the fermion necessarily has never been found 
before. As a summary, the respective coupling orders of 
the spectral properties of the ultrasoft mode are tabulated 
in Table. 1, together with those of the normal fermion and 
the antiplasmino obtained in the HTL approximation. 

The pole given by Eq. (16) gives rise to a new peak in 
the spectral function of the fermion as 



. Z 
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which is depicted in Fig. 3, where |p| is set to zero. Since 
the expression of £ is not available in the literature, we 
simply adopt £ = g 4 T In g -1 /(2ir) in the figure. 

3. Absence of Vertex Correction 

So far, we have considered the one-loop diagram. We 
need to check that the higher-order loops are suppressed 
by the coupling constant. This task would not be straight- 
forward because, 5m 2 ~ g 2 T 2 appears in the denominator, 
as seen in Eq. (13), which could make invalid the naive ex- 
pansion in terms of the coupling constant. The possible 
diagrams contributing in the leading order are ladder dia- 
grams shown in Fig. 5 because the fermion-boson pair of 
the propagators gives a contribution of order 1/g 2 , and the 
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Figure 4: The dispersion relation in the fermionic sector. The cou- 
pling constant is set to g = 0.1. The vertical axis is the energy ui, 
while the horizontal axis is the momentum |p|. The solid (blue) lines 
correspond to the normal fermion and the antiplasmino, while the 
bold solid (red) one to the ultrasoft mode. Notice that since we fo- 
cus on the fermion sector, the antiplasmino appears instead of the 
plasmino. The dotted lines denote the light cone. Since our anal- 
ysis on the ultrasoft mode is valid only in |p| <g g 2 T, the plot in 
|p| ^ 9 2 T ma y n °t have physical meaning. The residue of the an- 
tiplasmino becomes exponentially small for |p| 3> gT, so the plot of 
the antiplasmino does not represent physical excitation in |p| 2> gT, 
either. 



vertex gives g 2 . However, there is a special suppression 
mechanism in the present case with the scalar coupling. 
For example, let us evaluate the first diagram in Fig. 5, at 
small p. The self-energy is evaluated to be 
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Since there are four vertices and two pairs of the propa- 
gators whose momenta are almost the same, the formula 
would apparently yield the factor <? 4 x (5m 2 )~ 2 ~ g°, which 
is the order of unity. One can easily verify that this or- 
der would remain the same in any higher-loop diagram, 
so any ladder diagram seems to contribute in the same 
leading order as explained. However, the explicit evalu- 
ation of the numerators of the fermion propagators gives 
— f)f = \jk 2 — ftl 2 , which turns out to be of order g 2 ; 
this is because the internal particle is almost on-shell, k 2 , 
I 2 ~ g 2 . In the higher-order diagrams such as the second 
diagram in Fig. 5, the same suppression occurs. Thus, the 
contribution of the ladder diagrams giving a vertex correc- 
tion is absent in the leading order in the scalar coupling, 
and hence the one-loop diagram in Fig. 2 with the dressed 
propagators solely gives the leading-order contribution to 
the self-energy. We remark that a similar suppression oc- 
curs in the effective three-point- vertex at p ~ gT [16]. It 
should be emphasized that this suppression of the vertex 
corrections is, however, not the case in QED/QCD, where 
all the ladder diagrams contribute in the leading order and 




Figure 5: Some of the ladder class diagrams. The solid and dashed line are the dressed propagators of the fermion and the scalar boson, 
respectively. 



must be summed over [9, 13]. 

4. Summary and Concluding Remarks 

We have investigated the spectral properties of mass- 
less fermion coupled with a massless scalar boson in the 
ultrasoft momentum region (jp < g 2 T) at high tempera- 
ture. We have first indicated that a novel resummed per- 
turbation theory is needed beyond the conventional HTL 
approximation to get a sensible spectral properties in the 
ultrasoft region: For a consistent calculation in the rele- 
vant order of the coupling, we have shown that the use of 
the dressed propagators with the asymptotic masses and 
the width for the both fermion and boson is necessary and 
sufficient; although the vertex corrections due to ladder 
diagrams is apparently to contribute in the leading order, 
they turn out to be of higher order in the coupling and can 
be neglected. We have found that the resulting fermion 
propagator develops a novel pole yielding a peak in the 
spectral function in the ultrasoft region: Its pole position, 
width, and the residue arc obtained for the first time in the 
present work as a function of the asymptotic masses and 
damping rates; they are summarized and compared with 
those of the normal-fcrmion and (anti-)plasmino excita- 
tions obtained in the HTL approximation [2, 3] in Table. 1. 
Thus, we have established the existence of an ultrasoft 
fcrmionic mode at high temperature in the case of scalar 
coupling for the first time, although there were sugges- 
tions of the existence of such an ultrasoft fermionic mode 
[4, 5, 6, 9]. We remark that although the present analysis 
is given for a massless fermion, an ultrasoft fermion exci- 
tation should exist even for a massive fermion, if the bare 
mass is smaller than g 2 T, which is the smallest scale in 
our analysis. 

We should also emphasize that the existence of a new 
peak at the origin implies that the fermion spectral func- 
tion has a three-peak structure, although the central peak 
has only a small strength. The development of such a 
three-peak structure at intermediate temperature was sug- 
gested in the case where the boson is massive irrespective 
of the type of the boson [4, 5, 6], as mentioned in the In- 
troduction. In this sense, the three-peak structure of the 
fermion spectral function is a robust and persistent phe- 
nomenon to be seen at intermediate and high T, at least 
in Yukawa theories. 



Now let us discuss the physical origin of the ultrasoft 
mode for clarifying its possible universality or robustness 
at high temperature. As mentioned in the Introduction, 
an ultrasoft fermionic mode with a vanishing mass can 
appear as a phonino associated with the spontaneously 
broken supersymmetry at finite T [7, 8, 9]: In the mas- 
sive Wess-Zumino model, the supersymmetric cancellation 
of the fermion mass is essential to make the phonino ex- 
citation [8]. In fact, the massless fermionic mode is re- 
alized as a consequence of an exact cancellation of the 
self-energy with the finite bare mass at the vanishing ex- 
ternal momentum due to the spontaneous breaking of su- 
persymmetry. This is quite in contrast to the case of the 
Yukawa model dealt in the present work, where the ultra- 
soft fermionic modes do not appear when the fermion mass 
is large. This means that the ultrasoft fermionic mode 
like the phonino can appear even in a nonsupersymmet- 
ric model like the Yukawa model, and there the massless- 
ness of the fermion is an essential ingredient in realizing it. 
Nevertheless it is interesting that such an exotic fcrmionic 
mode can appear both in supersymmetric and nonsupcr- 
symmetric models with quite different mechanisms at high 
temperature. For getting further intuition into the possi- 
ble mechanism for realizing the ultrasoft fcrmionic mode 
in the Yukawa model, it is noted that Kitazawa et al. ar- 
gued that the level repulsion due to the Landau damping 
causes the three-peak structure in Ref. [4], though in the 
case where a massless fermion is coupled with a massive 
boson and thus the fermion-boson mass difference squared 
Sm 2 is nonzero. Indeed a finite 5m 2 plays an important 
role in the appearance of the ultrasoft mode even when 
the boson is massless because it ensures the smoothness of 
the self-energy at the origin: As discussed in Introduction, 
the real part of the self-energy vanishes at the origin from 
the symmetry for the massless fermion. The residue of the 
pole is proportional to (dm 2 ) 2 and vanishes if 5m 2 = 0. 
In the present work, the mass difference squared is a re- 
sult by taking into account the effects beyond the HTL 
approximation and found to be 5m 2 ~ g 2 T 2 . The mass 
difference in turn causes the smallness of the imaginary 
part: The imaginary part of the self-energy can originate 
from the boson emission, the Landau damping, and the 
imaginary part (C/,&) of the dressed propagators obtained 
beyond the HTL approximation. The former two contri- 
butions are found to be zero at the origin [4, 6] due to the 
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nonzero mass difference or kinematics at the one-ioop or- 
der. Therefore, the leading contribution of the imaginary 
part is solely given by the damping rates of the hard parti- 
cles of order g 4 T\ng~ 1 . As a result, the ultrasoft fermionic 
mode appears as the sharp peak at the origin with a small 
residue. 

Though the present analysis is limited to Yukawa the- 
ory with a scalar coupling, a similar analysis can be per- 
formed in other massless fermion-boson system including 
QED and QCD, although some complications arise [9, 13]: 
The vertex correction is not negligible in the gauge theo- 
ries in sharp contrast to the Yukawa theory, so a resum- 
mation of the ladder diagram is necessary. We also note 
that the damping rate of the hard fermion is of the or- 
der of g 2 T\ng~ 1 and "anomalously" large in QED and 
QCD [21]. Then the imaginary part of the fermion pole 
in the ultrasoft region should be larger than that in the 
Yukawa theory and the nature of the ultrasoft mode may 
be different from that in the present case [13]. In QCD, 
the self-coupling between the gluons should be also taken 
into account [9, 13]. All these matter will be reported in 
a separate paper [13]. 
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